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I. INTRODUCTION 

The characteristics of ultracold bosonic two-body scat- 
tering processes provide a key ingredient for modeling 
and understanding of the properties of degenerate atomic 
gases. The dominant s-wave interactions of the two-body 
collisions can be conveniently modeled by Huang-Fermi's 
pseudopotential [1-3]. Over the last decade, this pseu- 
dopotential has become an essential theoretical tool in 
particular for analytical studies of ultracold systems, e.g., 
see Busch et. al. for two atoms in harmonic traps, 
Olshanii for atomic collisions in cylindrical waveguides 
[5], where the so-called confinement-induced resonances 
(CIRs) was predicted and experimentally observed in [6,- 
[il, or in refs. (ol-ITTI which employed approaches based 
on the energy dependent scattering lengths, as{E). 

The theory of pseudopotentials has been developed fur- 
ther in order to describe scattering processes involving 
higher-partial waves, i.e., with relative angular momen- 
tum ^ > 0, i, [iiHioj. Higher-partial wave interac- 
tions represent an intriguing perspective in research on 
ultracold gases and are expected to provide novel many- 
body phenomena, such as unconventional superfluidity 
and superconductivity [2ll-[23|. So far, theoretical stud- 
ies with pseudopotentials for ^ > have mainly been 
used to describe fermionic systems, e.g., in quasi-one- 
dimensional pJI, [l5|, [2^ or quasi- two-dimensional [13, [l3 
geometries, which have been studied experimentally in 
[25| . For the case of bosons an analytical description of 
colHsions in a planar waveguide has been derived, where 
the cross-section possesses s- and d-wave resonances [17]. 
In harmonic waveguides higher-partial wave interactions 
were shown to have important impact on bosonic scat- 
tering [26] through the interference of d-wave resonances 
(DWRs) with the underlying continuum. A theoretical 
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treatment of d-wave scattering of bosons in the case of 
quasi- ID geometries is, however, still needed. 

In this work we derive an approximative analytical so- 
lution for bosonic collisions in a harmonic waveguide, 
which goes beyond mere s-wave scattering. Our starting- 
point is to assume that the transversal confinement fre- 
quency is weak enough in order to avoid a strong inter- 
ference between the free space d-wave shape resonance 
and the s-wave scattering. The latter permits to model 
the interatomic potential through a zero-range regular- 
ized d-wave pseudopotential suggested in p^,ll5|, which 
is parametrized by a single parameter, the energy depen- 
dent d-wave scattering length {ad{k)). The associated 
regularization operator contains exclusively derivatives. 
Taking into account the virtual excitations occurring dur- 
ing the scattering process we analytically describe the 
appearance of the DWR. An important aspect of our ap- 
proximative formula for the d-wave scattering amplitude 
is that one can map the total Hamiltonian onto an effec- 
tive one-dimensional Hamiltonian, where the interaction 
is modeled by a one-dimensional contact potential. The 
strength of the contact potential depends on the above- 
mentioned d-wave scattering amplitude, similar to the 
case of the s-wave CIR [5]. 

In order to verify the validity of the approximative an- 
alytical results we compare them with corresponding nu- 
merical calculations, where the interatomic interaction 
is modeled by a Lennard- Jones potential. A good agree- 
ment is obtained only for a sufficiently weak confinement. 
For strong confinement the pseudopotential does not cap- 
ture the s- and d-wave interference and therefore lacks the 
description of the numerically obtained Fano asymmetric 
lineshape of the transmission coefficient [27]. 

The sectioning of this paper is as follows. In Sec. II we 
review the recent developments of pseudopotential the- 
ory for higher partial wave interactions. In Sec. HI we 
describe and discuss the d-wave scattering process. Sec. 
IV contains a comparison of the analytical results and 
the numerical calculations and Sec. V provides a brief 
summary and conclusions. 
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II. PSEUDOPOTENTIALS FOR HIGHER 
PARTIAL WAVE INTERACTIONS 

The pseudopotential approximation to the interatomic 
interaction is an ubiquitous tool for the theoretical de- 
scription of ultracold gases: the atomic interactions are 
modeled in terms of a contact potential. Huang and Yang 
[2I, [3] derived a generalized pseudopotential for all par- 
tial waves, which represents a multipolar expansion over 
delta function contributions. The Huang- Yang approach 
(HY) was revised by Derevianko [12], who corrected an 
algebraic error for ^ > 0. However, other derivations have 
been developed for higher partial wave pseudopotentials, 
which differ from the HY-approach, by using a more rig- 
orous treatment of the derivatives of the delta-functions 
in the framework of distribution theory [3, [IH yielding 
the following representation: 



The above-mentioned pseudopotentials can be simpli- 
fied [15] in order to describe interactions of a single par- 
tial wave character. This simplification is based on the 
fact that in general the atomic interaction is dominated 
by the ^-partial wave scattering in the vicinity of the 
corresponding ^-partial wave resonance. Thus, one can 
focus on a single partial wave, which yields the following 
pseudopotential: 



T/p,,,(r)*(r) = 



167r2(-l)< 



geYim{d)S{r) 



Or 



(2£ + l)! 

r -V;„(a')*(r') 



2£+l ' 



(■4) 



r'=0 



Furthermore, the sum w.r.t. m can be simplified with 
the help of the following summation formula: 



X. , ^T / ^ (-1)M(2£+1)!!]' 5W(r) 



£=0 



dr 



47r(2^)!^! 



J r'=0 



(1) 



where = fi^af +^(A:)/(2/i), with af^\k) = 

— tan(5^(/c)//c^^+^ defined as the ^-th energy-dependent 
scattering length [lo|,[l^. Here, S^^\r) is the ^-th deriva- 
tive of the delta function and the integral over the solid 
angle fl^ acts as a projection operator of the total wave- 
function on each ^— state. For i = one obtains from 
Eq.([T]) the Huang- Fermi pseudopotential. However, the 
Eq.© for £ > differs by a factor {2£ + !)/{£ + 1) 
from Huang and Yang's result [3], whereas if we sub- 
stitute in Eq.([T]) the relations 5{r) = ^(r)/(47rr^) and 
= {—iy£\S{r)/r^ we obtain the pseudopotential 
of Derevianko [12]. 

Additionally, Stampfer and Wagner [16] introduce the 
following pseudopotential for higher partial waves: 



£=0 m=- 



47r(-l) 



-g£Yim{d)S{r) 



X 



{2£)\\ 



r'=0 



(2) 



where the operator Yg^(9), introduced by Maxwell [28|, 
is obtained from the spherical harmonic polynomial 
r^Yimi^)^ by replacing the Cartesian coordinates Xk with 
the partial derivatives dx,^. Note that Eq.(|2]) is equal to 
Eq.(P) as it was proven in by using the relation 



47r£! 



{2£+l)\\ 



Ye„,{d)5{r). 



(3) 
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[e/2] 



2i- 

m=-t /c=U 

(5) 

where [£/2] is the largest, smaller or equal integer of 
£/2 and cu = [(-1)^(2^ + 1)!] / [2^i^!(^ - - 2A:)!] 
are the coefficients of the Legendre polynomials Pi{x) = 

Eqs.dlJ and (|5|) constitute pseudopotentials for any i- 
partial wave. These pseudopotentials are expressed in 
terms of differential operators which are more convenient 
in analytical derivations than the pseudopotential with 
the projection operator (see Eqs.([T]) and (|2])). 



III. QUASI-ID SCATTERING WITH S- AND 
D-WAVE INTERACTIONS 

A. s- and d-wave pseudopotentials 

In the following we consider ultracold collisions of iden- 
tical bosons in the presence of an external harmonic 
confinement. The scattering process is taking place in 
the low energy regime, ka± <C 1, where k is the rela- 
tive momentum in the unconfined degree of freedom and 
= oscillator length defined via the trap 

frequency uj± and the reduced mass /i of the colliding 
bosonic pair. Furthermore, in order to go beyond s-wave 
physics we assume that the momentum k is very small 
but not equal to zero. Consequently, as it was shown 
in [26] except for the well known s-wave CIR [5], an ad- 
ditional resonance emerges in this system, the so-called 
DWR, which is based on the interference of a free space 
d-wave shape resonance with s-wave scattering. The co- 
existence of s- and d-wave resonances implies interactions 
between the bosons of s- and d-wave symmetries (see 
Fig. 1(c) in Ref.[26]). 

The main goal of this paper is to present an analyti- 
cal model for the DWR effect in a harmonic waveguide. 
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However, the modeling of s- and d-wave interactions with 
pseudopotentials in systems of non-spherical symmetry, 
as the above-mentioned one, is not a trivial task. As 
the symmetry of the system is cylindrical the angular 
momentum is not conserved and consequently, the appli- 
cation of pseudopotentials, Eqs. ^ and (j2j), leads to an 
explicit integration over all ^-partial waves. 

Though, as it was proposed in [15] one can avoid 
such difficulties by considering a single -^-partial wave in 
the vicinity of the ^-wave resonance, where contributions 
from other partial waves are sufficiently suppressed. The 
latter means that in the vicinity of the s-wave resonance, 
we assume that the interaction can be modeled by the 
following pseudopotential with the help of Eqs. ^ and 
®: 



^^fiulp^^VpsA^), (8) 

where r = ri — r2 is the relative coordinate of the two 
bosons, and Vps/{r) is the pseudopotential which models 
the arbitrary ^-wave inter-atomic interaction. 

The symmetry of the transversal potential, ^/icj^p^, 
implies that the scattering solutions should be expanded 
in an axially symmetric basis. Such a basis set are the 
eigenstates of the two-dimensional (2D) harmonic oscil- 
lator which satisfy the Schrodinger equation 



27Th^as{k) d 
VpsA^) = S{r) — [v 



(6) 



which is the Huang-Fermi pseudopotential, with the s- 
wave scattering length as{k) being energy dependent. 
Note that Eq.® coincides with the pseudopotential pro- 
posed in [11]. 

Equivalent ly, in the vicinity of the d-wave resonance, 
we assume that the interaction can be modeled by the 
following pseudopotential [13, [HI ' 



8/i 



E 

i,j,k,l 



(7) 

where Dijki = SikSji — \5ij5ki with the indices z, j, /c, / = 
1, . . . , 3 and a^(A:) = — tanj£=2(^)/^^ is the energy depen- 
dent d-wave scattering length. The bidirectional opera- 
tor, ^ {d'^.^ . ) {{d'^^^^ )^) denotes the differential operator 
in Cartesian coordinates, that acts to the left (right) of 
the pseudopotential VpsA^)- 

We emphasize that application of the expressions (|6]) 
and ([7j) simultaneously is inconsistent and incompatible 
since the pseudopotentials lps,s(r) and VpsA^) valid 
for the cases where the wavefunction exhibits only 1/r- 
or 1/r^- singularities for small r, respectively. 



Hamiltonian, scattering states and virtual 
excitations 



We assume two-body collisions of identical bosons in 
a waveguide with their transversal degrees of freedom 
being confined by a harmonic oscillator potential. The 
harmonic confinement permits the separation of the cen- 
ter of mass and relative motion yielding a single-particle 
Hamiltonian with a scatterer fixed at the origin for the 
relative degree of freedom. 



2 II \dp^ 



i_a 

pdp 



1 



L 2 2 

-puj^p 



(9) 



where ^n.m^ {p-, 0) are the eigenfunctions of the 2D har- 
monic oscillator, is the azimuthal quantum num- 
ber and En^m^ = hw±{n + \mz\ + 1) is the transver- 
sal energy spectrum of the 2D harmonic oscillator, with 
n = 0, 2, 4, . . . being the principal quantum number. 

The spherical symmetry imprinted by the pseudopo- 
tential, Vps/{t)^ leads to a separation of the azimuthal 
(j) angle from the z, p coordinates in Eq.(|9]). The latter 
means that during the collision no virtual excitations will 
occur with respect to the quantum number. However, 
virtual excitations will emerge with respect to the princi- 
pal quantum number n, which will be taken into account 
in the calculation given below. 

We focus on the single mode regime, where the col- 
lision takes place in the ground state of the transversal 
confinement with angular momentum = and repre- 
sents the only "open channel" . The latter means that the 
total energy of the pair collision is limited between the 
ground and the ffist axially symmetric transversal state 



n=z,mz 



(10) 



Furthermore, the above-mentioned independence on 
the angle <p results in the following simplified Schrodinger 
equation: 



^2 / 5^ 



^ d\ 1 2 2 T. / 



^{z,p) 



E^{z,p), 
(11) 



where E = E\\ -\- En=o^rnz=o is the total colliding en- 
ergy, which is defined as the sum of the longitudinal ki- 

netic energy E\\ ~ energy of the transverse 
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ground state En=o^m^=o = huj±. Due to the constraint of 
Eq. (fTQ|) the momentum k is limited according to the rela- 



tion k < 



2/i 



=o). Here, ^(z,/?) 



denotes the full 2D axially symmetric solution, which can 
be expanded in the 3n,o{p^ 0) = ^n{p) basis 



(12) 



n=0 



where we sum over all even n due to the bosonic symme- 
try. 

Substituting Eq. (p!2]) into Eq. (pT]) , multiplying by 
27rp(?*(p), integrating and using Eq.(j9]) yields the fol- 
lowing equation for the functions Cn{z): 



E 



f 

Jo 



(13) 



By solving Eq. (p!3|) for n = we obtain the scattering 
wavefunction Co{z) in the open channel which has the 
following asymptotic form: 



Vd = 



with p given by the relation p = y^x^hT^. 

Then the expression for the wavefunction reads 



^{z,p) 



48al 



where Ai(z,p) is defined as follows: 



A / , ^ 2n + l 



^0+ 

(17) 



ridAi{z,p), (18) 



(19) 

Let us now to proceed with the explicit definition of 
the regular part rjd of the wavefunction, being a solution 
of Eq.([T7j). The action of the operator [d^ + S^) on ^(r) 
for X, ^ results in the following relation: 



Co(^) =cos(fcz) + /e,£e^^l^l, for|z|^oo, (14) 

where fe/ is the one-dimensional scattering amplitude of 
the ^-wave interaction which describe the even scattered 
waves. 

For n > 0, Eq. (p!3|) provides us with the solution for the 
virtual excitations associated with the closed channels, 
which decay exponentially according to the relation 



Cn{z) = An/ 



2\z\ 



(15) 



where the coefficients An/ refer to the ^-wave interaction 
and denote the transition amplitudes from the transverse 
ground state to the n-th excited state. 

In the case of the d-wave pseudopotential we substitute 
Eq. ([T4|) into Eq.([l3l) for n = and Eq.l^ into Ea. ffn]) 
for n > and integrate over the z variable , respectively, 
in the interval [— e, e] with e ^ 0. Finally, we obtain the 
following relations for the d-wave scattering and transi- 
tion amplitudes: 



.al{k)^ 

Je,d ^ 



24kal 



Vd ; ^n,d 



An,_d — 



(^di^)V^ 2n + l 



48al 



-Vd 



n / fca_ 

2 V 2 



2 

(16) 



where rjd is the regularized part of the wavefunction ^ in 
the limit r ^ 0, 



{dl + dl)^{z,p) 



■ ^dWV^ ^ Jk\z\ 



24a\ 



r]d^2{z,p = 0), 
(20) 



where A2(z, p = 0) reads: 



A2(z,p = 0) = ^ 



21^1 



(21) 

Here, we should note that the A-series at p = are 
not converging uniformly as z ^ 0. The latter means 
that the partial derivative w.r.t. z-coordinate cannot be 
interchanged with the sum over all even n. Thus, firstly 
one has to perform the summation and then differentiate 
with respect to the z coordinate. However, with the help 
of Euler's summation formula [29j one can express the 
A-series as an expansion with respect to the z variable, 
where the singular part of the A-series will be written in 
closed form yielding the following expressions: 



2a\ 



2\z\ 



|^|3 



2\z\ 



a± 



A2 z, p = = —-^ + ^ + ^ + ^1 



(22) 
(23) 
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where qi ^ -2.29225, ^2 = 0.83333, qs ~ -0.154912 
and bi ^ —3.5336 are constants calculated via Euler's 
summation formula. 

Consequently, by substituting Eqs. (p!8|) - (|23|) into 
Eq. (fT7|) we obtain r]d^ which is now free of divergences. 
This in turn yields the following expression for the one- 
dimensional scattering amplitude: 



fe,d 



ka± \ ( _ 



10a^(/c) 



(24) 



where the constant 62 ^ 2.386445 is defined by the re- 
lation 62 = — "Y — 4^3. The resonance condition for f^^d 
then reads 



ad{k) 



\/i06^. 



(25) 



Eq. ()24|) shows that the d-wave scattering in the pres- 
ence of the waveguide becomes resonant at an off- 
resonant value of the free space ad{k) scattering length, 
due to the significant contributions of the virtual exci- 
tations. The latter holds equally for the case of s-wave 
CIR. Therefore, one can conclude that in general the ap- 
pearance of resonances in systems with transversal con- 
finement are generated by the virtual excitations in the 
transverse modes. 

By applying the s-wave pseudopotential, given in 
Eq.(j6]), and performing the derivation as mentioned 
above, we will obtain the same relation for the one- 
dimensional scattering amplitude as in [5] with the dif- 
ference that in this case the s-wave scattering length is 
energy dependent [9,]: 



/e,s — 



1 



1 + i 



^^^^ + 



(26) 



The constant ci in Eq. ([26|) is equal to 1.4603. 



IV. RESULTS AND DISCUSSION 

In this section we compare the analytical results of the 
previous subsection for the s- and d-wave pseudopoten- 
tials with the corresponding numerical simulations. First 
we evaluate the d- and s-wave energy dependent scatter- 
ing lengths by solving numerically a model of Cs atoms 
in free space interacting via a Lennard- Jones potential, 
V{r) = — Note that Cs atoms are ideal canditate 
for studing scattering phenomena, since they provide a 
rich spectrum of resonances [30|, . 

The dispersion coefficient Cq has been taken from 
[33I and C12 is a free parameter, which controls the 
values for the s- and d-wave scattering lengths. This 
yields a parametrization of the one-dimensional s- and 



d-wave scattering amplitudes, Eqs. pH]) and (|24|) . respec- 
tively, in terms of C12. 

In order to numerically solve the Cs-Cs collisions in 
the harmonic waveguide we employ the units mcs/2 = 
h = ujQ = 1, where mcs is the mass of the Cs atom and 
cjo = 27r X 10 MHz. The longitudinal energy is set to e\\ = 
2 X 10~^ and the transversal energy is varied within the 
interval 10~^ < < 8 x 10~^, corresponding to a range 
0.27r kHz < uj± < 167r kHz for the waveguide confinement 
frequency. We thereby focus on the low energy regime, 
characterized by ka± <C 1. In the following we calculate 
the transmission coefficient given by: 



Te = \l + fe,ef, 



(27) 



where fe/ is the one-dimensional scattering amplitude 
for the ^-wave. 

In Fig.[T]we show the transmission coefficient as a func- 
tion of C12 with the transversal frequency uj± = 10~^. 
Tn in Fig. [T]^a) denotes the numerically calculated trans- 
mission coefficient, as it was shown in [26, 34, 35]. We 
observe the appearance of two minima (T^ = 0) due to 
resonances. The minimum on the r.h.s. of Fig.[TJa) refers 
to the s-wave CIR and the minimum on the l.h.s. to the 
DWR. Due to the centrifugal barrier the width of the 
DWR is much narrower than that of the s-wave CIR. 
Additionally, we observe that both resonances yield a 
Lorentzian-like lineshape in the transmission spectrum. 
This occurs due to the fact that the weak confinement 
suppresses the interference effects between the s-wave and 
d-wave interactions. 

In Fig.fTJb) the analytically calculated transmission co- 
efficient Tg, as given by Eqs. (|26l) and (|27|) . is shown which 
refers to s-wave interactions only. Tg describes very accu- 
rately the s-wave CIR (see Fig. [Ha)-(b)), whereas there 
is no trace of the DWR, since contributions from higher 
partial waves are not included in Eq. (|26l) . In Fig. [TJc) 
the analytic expression of the transmission coefficient Td 
is depicted according to Eqs. (|24|) and ([27]) . We observe 
that the DWR is described qualitatively accurate enough 
(see also insets of Figsl2fa)-(b)), despite the fact that in 
the analytical calculations we have neglected the contri- 
butions of s-wave scattering, since the transversal con- 
finement is very weak. We note that the results shown 
in Figs.[Hb) and[Hc) being obtained analytically for the 
pure s- and d-wave scattering cases lead to an accurate 
description of the scattering processes in the presence of 
a weak external harmonic potential. 

Let us now determine the regime of validity of the ap- 
proximative analytical results for the d-wave pseudopo- 
tential scattering by comparing it with the corresponding 
numerical simulations. In Fig. [2] we show the transmis- 
sion coefficient T calculated numerically (Fig. [2fa)) and 
analytically (Fig. [2fb)) for several transversal confine- 
ment frequencies uj±. In Fig. [2fa) for a strong confine- 
ment, e.g. uj± = 8 X 10~^, we first observe that the 
numerically calculated position of the resonance deviates 
substantially from the analytical one. However, as uj± de- 
creases the numerical results for the position of the DWR 
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FIG. 1: (a) Numerically calculated transmission coefficient 
Tn, which shows the DWR on the l.h.s and the s-wave CIR on 
the r.h.s.. The analytically calculated (b) s-wave transmission 
coefficient, Ts, and (c) d-wave transmission coefficient, T^. 
The transversal confinement is LjJ± = 10"^ 



converges to the corresponding analytical one. The ori- 
gin of this behavior is the fact that we considered only 
d-wave interactions in the analytical calculations and ne- 
glected contributions from the s-wave interactions. This 
approximation is eligible in the regime of a weak confine- 
ment, e.g. uj± ^ 4 X 10~^, where the s- and d-wave free 
space resonances are weakly coupled, and in the vicin- 
ity of the DWR the d-wave interactions are dominant 
over the s-wave. In the regime of strong confinement the 
transmission spectrum T, in the vicinity of the DWR, 
exhibits a strong asymmetric Fano-lineshape, with Tn 
changing abruptly from zero to one. The latter is an 
interference effect of the strong s- and d-wave interac- 
tions in the confinement and cannot be described by this 
analytical pure d-wave model. However, this asymmet- 
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FIG. 2: (color online) (a) (left panel) The numerically calcu- 
lated transmission coefficient Tn versus C12 parameter for sev- 
eral confinement frequencies uj± , (right panel) high-resolution 
figure of the DWR for a;^ = 10~^ where the inwards-pointing 
arrows denote the width of the resonance. (b)(left panel) The 
analytically calculated transmission coefficient Td as a func- 
tion of the C12 parameter calculated with the d-wave pseu- 
dopotential for the corresponding values of uj±, (right panel) 
high-resolution graph of the DWR for ul)± = 10~^. The yel- 
low dashed line in the right panels indicates the minimum of 
the numerical calculated transmission coefficient (Tn = 0) in 
comparison with the corresponding analytical result. 



ric profile gets suppressed with decreasing uj_\_ yielding a 
Lorentzian-like lineshape. 

The right panels of Fig. [2] show a high-resolution graph 
of the DWR calculated analytically and numerically for 
uj±_ = 10~^, where the dashed vertical line in both plots 
indicates the position of the minimum of the numerically 
calculated transmission coefficient T^. We observe that 
the d-wave pseudopotential approximation qualitatively 
describes the numerical simulations in the regime of weak 
confinement. However, in this small scale of the C12 
values the analytical results show that the width of the 
DWR is by two orders wider than the numerical predic- 
tions, as well as there is a slight difference in the position 
of the DWR. These quantitative deviations again occur 
due to s-wave interactions, which act as corrections in the 
transmission coefficient Tn despite of the weak coupling 
of the s-wave and d-wave interactions. 



V. BRIEF SUMMARY 

We have derived and analyzed analytical expressions 
for resonant d-wave scattering in a harmonic waveguide, 
in the framework of pseudopotential theory. The inter- 
atomic interaction has been modeled by a pure d-wave 
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pseudopotential introduced in [13, [HI . We observe that 
the analytical pseudopotential approach provides results 
which are in qualitative agreement with the correspond- 
ing numerical results in the regime of weak confinement 
due to the very small interference effects of s- and d-wave 
inter act ions. The appearance of the DWR in a waveguide 
is strongly affected by the virtual excitations over all the 
closed excited modes of the transversal confinement. For 
a strong transversal confinement the analytically and nu- 
merically calculated lineshape and position of the DWR 
deviate from each other. The interference effect of the 
s- and d-wave free space resonances in the confinement 
yields an asymmetric Fano profile which cannot be de- 
rived by a single partial wave model. Our analytical re- 
sult Eq. (|24|) for weak confinement can be used for the 
analysis of resonant scattering processes in confined ge- 
ometries beyond s-wave physics. Finally, we remark that 
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